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Abstract 

This paper introduces an interpolation framework for the weighted-?^ model reduction problem. We obtain a new representa- 
tion of the weighted-?^ norm of SISO systems that provides new interpolatory first order necessary conditions for an optimal 
reduced-order model. The H2 norm representation also provides an error expression that motivates a new weighted-?^ model 
reduction algorithm. Several numerical examples illustrate the effectiveness of the proposed approach. 
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1 Introduction 

Consider a single input/single output (SISO) linear dy- 
namical system with a realization 

Ex(t) = Ax(i) + bu(t), y(t)=c T x(i) (1) 

for E, A e R nxn and b,c E W\ x(i) E R n , u{t) E E, 
y{t) E M, are respectively the state, input, and out- 
put of the system. The transfer function of this system 
is G(s) — c T (sE — A) _1 b. Following common usage, 
the underlying system will also be denoted by G. For 
many examples, the state-space dimension n is quite 
large, leading to untenable demands on computational 
resources. Model reduction attempts to address this by 
finding a reduced-order system of the form, 

E r x r (t) = A r x r (i) +b r u(t), y r {t) = c^x r (t) (2) 

with G r (s) = cJ{sE r - A,.)- 1 ^ for E r , A r E R rxr 
and b r , c r E W with r <C n such that y r (t) w y(t) 
over a large class of inputs u(t). G r is a low order, yet 
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high fidelity, approximation to G. We construct G r 
via state-space projection: Two matrices ("reduction 
bases") V r , W r eR" xr are chosen. Then, system dy- 
namics are approximated by x(t) w V r x r (i) and forcing 
a Petrov-Galerkin conditon ("orthogonal residuals") 
Wj(EV r x r (t) - AV r x r (i) -bu(t)) = 0, together with 
the output equation y r (t) = c T V r x r (t) to produce 

E r = W r T EV r , b r = W^b, 
A r = WjAV r , and = c T V r . 

See [2,3] for more information. 

1.1 Model Reduction by Interpolation 

The reduction bases, V r and W r , used in (3) will be 
chosen to force interpolation: G r (s) will interpolate G(s) 
(possibly together with higher order derivatives) at se- 
lected interpolation points. This approach to rational 
interpolation has been considered in [20,21,5,8,7,3] and 
depends on the following result. 

Theorem 1 Given two sets of interpolation points 
{°~k}k=i and {£fcK = i; that are each closed under con- 
jugation, and a dynamical system G as in (1), consider 
matrices V r and W r such that 

RangeCVr) = span {(iTiE — A) _1 b| 

for i — 1, . . . , r. (4) 

Range(W r ) = span {((iE — A) c} 
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Then, V r and W r can be chosen to be real; G r (s) = 
c^(sE r - A,.)- 1 ^ defined by (2)-(3) is a real dynami- 
cal system that satisfies G{a k ) = G r (a k ) and G(( k ) = 
G r ((k) for k — 1, . . . , r; and, ifuj = Q for some j, then 
G'(o-j) = G' r (<Tj), as well where G' denotes the derivative 
ofG(s) with respect to s. 

Theorem 1 can be generalized to higher-order derivative 
interpolation as well, see [20,21,5,8,7,3]. The subspaces 
of Theorem 1 are rational Krylov subspaces and so, in- 
terpolator model reduction methods for SISO systems 
are sometimes referred to as rational Krylov methods. 



1.2 Weighted Model Reduction 

The Woo norm of a stable linear system associ- 
ated with a transfer function, G(s), is defined as 
||G|| Woo = max\G(iuj)\ . The "H 2 norm of G is defined 

1/2 



as \\G\\ n2 



\G(iLo)fdio 



The vector 



spaces of meromorphic functions that are analytic in 
the right halfplane, having cither bounded "Hoo norm or 
bounded H2 norm will be denoted simply as Hoc or "H 2 , 
respectively. Let W € "Hoo be given. The (W-)weighted 
%2 norm is defined as ||G||« 2 ( W) = \\G ■ W\\u 2 - 

We are interested in finding a reduced-order model G r 
that minimizes a M^-weighted TL2 norm, i.e., that solves 

\\G — G r \\n 2 (w) = min \\G - G r \\u 2 (w) (5) 

dim(G r )=r 

The introduction of W(s) allows one to penalize the error 
in certain frequency ranges more heavily than in others. 

An illustrative example: controller reduction 

Consider a linear dynamical system, P (the plant) to- 
gether with an associated stabilizing controller, G, that 
is connected to P in a feedback loop. Many control 
design methodologies, such as LQG and % 00 methods, 
lead ultimately to controllers whose order is gcnerically 
as high as the order of the plant, see [17,22] and refer- 
ences therein. Thus, high-order plants generally lead to 
high-order controllers. However, high-order controllers 
are usually undesirable in real-time applications due to 
complex hardware, degraded accuracy, and degraded 
computational speed. Thus, one prefers to use a reduced 
controller G r to replace G. Requiring G r to be a good ap- 
proximation to G is often not enough in terms of closed- 
loop performance; plant dynamics need to be taken 
into account during the reduction process. This may be 
achieved through frequency weighting: Given a stabiliz- 
ing controller G, if G has the same number of unstable 
poles as G r and if || [G - G r ]P[I + PG}^ 1 \\ n ^ < 1, 
then G r will also be a stabilizing controller [1,22]. Hence 
the controller reduction problem may be formulated 



as finding a reduced-order controller G r that mini- 
mizes or reduces the weighted error || (G — G r )W\\ n 
with W(s) := P(s)(I + P^Gis))- 1 ; i.e., controller 
reduction becomes an application of weighted model 
reduction. This approach has been considered in 
[17,1,14,10,6,19,12,18,16] and references therein, leading 
to variants of frequency-weighted balanced truncation. 
Conversely, the methods in [11] and [15] are tailored in- 
stead towards minimizing a weighted-?^ error as in (5). 

2 Weighted-?^ model reduction 

The methods proposed in [11] and [15] for approaching 
(5) require solving a sequence of large-scale Lyapunov 
or Riccati equations; they rapidly become computation- 
ally intractable as the system order, n, increases. We 
will approach (5) within an interpolatory model reduc- 
tion framework requiring only the solution of (generally 
sparse) linear systems and no need for dense matrix com- 
putations or solution of large-scale Lyapunov or Riccati 
equations. Interpolatory approaches can be effectively 
applied even when n reaches the tens of thousands. 

2.1 A representation of the weighted-V.2 norm 

Given transfer functions G, H e H2, and W G Hca, 
define the weighted-?^ inner product as 

(G, H) H2(W) = -L f°° G{iu)W{iu)W{iu)H{ico)du 
1 f x 

= — — / G{-lU})W{-lU})W{lU))H{%L0)du, 

2tt J -00 



so that \\G\\h 2 {w) = \J{G, G) U2 

(w)- The following 
lemma gives a compact expression for the weighted-?^ 
inner product based on the poles and residues of G(s), 
H(s) and W(s). By res[M(s), tt], we denote the residue 
of M(s) at tt e C. ' 

Lemma 2 Suppose G, H € H2 have poles denoted re- 
spectively as {Ai, . . . , A„} and{fii, . . . , /i m } 7 and suppose 
W € Hoc has poles denoted as {71, . . . ,7 P }. Assume that 
H(s) and W(s) have no common poles, and the poles of 
W(s) are simple. Then 

m 

fe=i 

+ EG(- 74 )^(-7,)#(7*) ' res[W(*),7i] 



Define X k = res[G(-s)W(-s)W(s)H(s), 
• If is a simple pole of H(s), then 

Xk=G{-Hk)W{-ix k )W{n k ) ■ r&s[H{s),n k }. 
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If /ik is a double pole ofH(s), then 

Xk =G(-fi k )W(-n k )W(^ k ) ■ res[H(s),fi k ] 

-h-M ■ ^ s [G(s)W(s)W(-s)]\ s= _ n 

where h-2{n k ) = lim ( s — fJ-k) 2 H(s). 



Proof: T(s) = G(-s)W(-s)W(s)H(s) has poles at 

{-Ai, . . . - A„} U {±71, . . . , ±7 P } U {(n, . . . , fi m }. 

For any R > 0, define a semicircular contour in the 
left halfplane: r^j = {z \ z = ius with ui <G [-R, R] } U 
{z \z = Re 10 with 6 e [f , ^f] } . For i? large enough, 
the region bounded by contains {71,..., 7 P } U 
{fj,i, . . . , Hm}i constituting all the poles of W(s)H(s), 
and hence all the stable poles of G(— s)W(— s) W(s)H (s) . 
Then, the Residue Theorem yields 

(G, H) H2(W) = ±- pGt-H^t-H^H^H^ 

J — 00 



= lim 



G(-s)W(-s)W(s)H(s) ds 



= res[G(-s)W(~s)W(s)H(s),u k ] 

v 

+ Y / res{G(-s)W{-s)W{s)H(s),-f i ]. 



This leads to the first assertion. If fi k is a simple pole for 
H(s), then 

r es[G(-s)W(-s)W(s)H(s), » k ] 

= lim [(s - (i k )G(-s)W(-s)W(s)H(s)} 

= G{-^ k )W{-n k )W{n k ) lim (a - n k )H(a). 

s->ftk 

Similarly, if /i k is a double pole for H(s), then it is also 
a double pole for G{-s)W{-s)W(s)H(s) and 



\G\\n 2 ( W ) 



]T G(-A fc )W(-A fc )W(A fc ) • res[G(s), X k ] 
fc=i 

p 

+ ^G(- 7fc )^(- 7fc )G( 7fc ) • res[^( S ), 7fe ]. (6) 



fe=i 

This new formula (6) for the weighted-?^ norm contains 
as a special case (with W(s) = 1), a similar expression 
for the (unweighted) H2 norm introduced in [9] . 

Suppose W € "Hoo has simple poles at {71, . . . , 7 p } and 
define a linear mapping $ ■ R2 —> R2 by 

m(s) = G(s)W(s)W(-s)+ir G{- lk )W{- lk ) ™Wifi)M 

k=l 8 + 7fc 

(7) 

Notice that G(s)W(s)W(— s) has simple poles at 
-7i,-72,---,-7 P , and 



res[G(s)W(s)W(-s),- lk } 

= lim (s + lk )G(s)W(s)W(-s) 

= G(- 7 fc)W(-7fc) I™ (s + l k )W(-s) 

s->-7fc 

= -G(- 7fc )W(-7fc) lim (s - l k )W{s) 

= -G(- 7 fc)W(-7fc)-res[W(a),7fc]. 

Thus g'[G](s) has poles only in the left half plane and 
indeed £ : "H 2 ^ "H 2 . 

Corollary 4 Suppose G and W are stable with poles 
{Ai, . . . , A„} and {71, . . . , j p }, respectively. Choose /1 ar- 
bitrarily in the left half plane distinct from these points. 

ThenforF(s) = d[G](s), (g, = F(-(j,) and 

Proof: By Lemma 2, 
G, — — 



t 1 / H 2 (W) 



G(-u)W(-u)W(u) 



+ Y,G{- lk )W{- lk ) 



res[VI^(s),7fc] 
7fc - A* 



= F(—fi), and 



res[G(-s)W(-s)W(«)IT(*), Mfe ] 

= lim A[ (s _ iUfe) 2 G( _ s)W( _ s)W(s) ^ (s)] 
s^pic as 

= lim G(- S )W(-s)W(s)^ [(* - ^k) 2 H{s)} 

+ lim (s - ^ k ) 2 H{s)^- [G{-s)W{-s)W{s)] 
s^Mfc as 

= G(-^ k )W{-,ik)W(pLk) ■ res[H(a), n k ] 
-h- 2 (Hk) ■ ^ [G(s)W(s)W(-s)]\ s= _ llk □ 

Corollary 3 If G(s) and W(s) in Lemma 2 each have 
simple poles, then 



G, 



- [G(s)W(s)W(-s)]\ s= _^ 



( s A 1 ) 2 / H 2 (W) 
+ ±G(- Jk )W { -, k ) r ^ V ^ = -F'(-,).n 



2.2 Weighted-H2 optimality conditions 

Consider the problem of finding a reduced order system, 
G r , that solves (5). The feasible set for (5) is noncon- 
vex, so finding a true (global) minimizer is generally in- 
tractable. Nonetheless, we are able to obtain descriptive 
necessary conditions for G r to satisfy (5). 
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Theorem 5 If G r has simple poles, {Ai, . . . , A r } ; and 
solves (5), then G r must satisfy: for k = 1, . . . , r, 

F r (-X k ) = F(-A fc ) and i?(-A fc ) = F'(-A fe ) (8) 
w/iere F = $[G] and F r = ${G r ] are defined from (7). 



Proof: Suppose by way of contradiction that, for some 
li e {Ai, A r }, (G-G r , ^) =a a ^0. By 

hypothesis, G r can be represented as G r (s) — J2l=i 



A, 



and for some index k, fi = A&. Define ?9 = ar g( a o) and 
with e > 0, define 



G^(s) = 

Then 

\\Gr-Gi e) \\u 2 (w) = 



s — fl 



E 



Ai 



S — fl 



<\\w\\ 



H 2 (W) 



so that ||G r (s) - Gr- e) (s)||« 2 (w) = O(e) as e 0. Since 
G r solves (5), 

||G-G r ||^ 2(M0 < ||G - G^||^ 2(lv) 
<\\(G - G r ) + (G r - G^)\\ 2 n2(w) 



<\\G-G r U2 



H 2 (W) 

+ \\G r -G^\\l 



2Re ( G - G r , G r -G< £) 

\ /«2(W) 



2(W)" 

Thus, 

< 2 Re (G - G r , G r - G< e >) +||G r -G< e >|& a(w0 . 

This implies first that < — e|«o| + C( £2 ): which then 
leads to a contradiction, olq = 0. 

To show the next assertion, suppose that for some p, € 
{Ai, . . . , A r }, ( G — G 



= ai ^ 0. Then 

IU 2 (W) 

for some fc, — X k and we define $i = arg(<p fe • ai). For 
e > sufficiently small, define 



G^( S ) = - 

As e — > 0, we have 
j|G r — Gl e '||« 2 (w) : 



E 



¥>i 



s - A, 



(s - n) 2 -ee~ 1 ^ 



= 0(s) 



H 2 (W) 



Following a similar argument as before, we find that 
< — e\0k • cki | + 0(e 2 ) as e — > 0, which leads to a con- 
tradiction, cti = 0. □ 



The interpolation conditions described in (8) give first 
order necessary conditions for G r to solve the optimal 
weighted- "H 2 model reduction problem (5). Note that for 
W(s) = 1 , one obtains F(s) = G(s) and F r (s) — G r (s); 
thus (8) contains the interpolatory "H 2 optimality condi- 
tions of [9] for the unweighted problem as a special case. 
Unfortunately, there does not appear to be a straight- 
forward generalization of the corresponding computa- 
tional approach that was described in [9] for the optimal 
(unweighted) H2 model reduction problem. Instead, we 
consider a different systematic approach to this problem 
motivated by an expression for the weighted- "H 2 error. 

2.3 A weighted-V.2 error expression 

The weighted- H2 norm expression in Corollary 3 leads 
immediately to an expression for the weighted-?^ error 
that forms the basis for our computational approach. 

Corollary 6 Suppose that G, G r and W are stable 

with simple poles {Ai}™ =1 , \ Xj> , and {7fc}^ =1; re- 

L J j=i 

spectively, and that there are no common poles. Define 
residues: &:= res[G(s), Xi}; 4>j:= res[G r (s), Xj]; and 
^pk-~ res[W(s),jk]- The weighted-V.2 error is given by 



G Gr 



H 2 (W) 



= £)((?(-Ai) - G r (-Ai))W(-Ai)W(Ai) • 4>i 



+ X)(G r (-A j ) - Gi-X^Wi-X^wCXj) ■ fa (9) 

p 

+ X)(G(-7fc) - G r {- lk ))W{- lk ){G{ lk ) - G r ( 7 *)) ' V>fe 
k=i 

One may recover the (unweighted) H2 error expression 
of [9] as a special case by taking W(s) = 1. Notice that 
the weighted error depends on the mismatch of G and 
G T at the reflected full system poles {— Ai}, reflected 
reduced poles {—Xj}, and reflected weight poles {— 7fc}- 

2.4 An algorithm for the weighted-V.2 model reduction 
problem: W-IRKA 

In order to reduce the weighted error, one may eliminate 
some terms in the error expression, by forcing interpola- 
tion at selected (mirrored) poles. Since r is required to 
be much smaller than n, there is not enough degrees of 
freedom to force interpolation at all the terms in the first 
and third components of the weighted- H2 error. How- 
ever, the second term, i.e. the mismatch at Xj, can be 

completely eliminated by enforcing G(— Xj) = G r (—Xj) 
for j — 1, . . . , r. Hence, as in the unweighted H2 prob- 
lem, the mirror images of the reduced-order poles play a 
crucial role. This motivates an algorithm with iterative 
rational Krylov steps to enforce the desired interpola- 
tion property as outlined in Algorithm I below. How- 
ever, a crucial difference from the unweighted H2 prob- 
lem is that we will not enforce interpolation of G'(s) at 
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Algorithm 1. Weighted Iterative Rational Krylov 
Algorithm (W-IRKA) 

Given G(s) = c T (sE - J 4)~ 1 b and VK(s) = c^(sE TO - 
A u ,) _1 bi 1) , reduction order r = + zu with > 0, 
let denote the v dominant poles of G and 

{ifejfcLi the dominant poles of W . 

(1) Make an initial interpolation point selection: 

Ci = -A, for i = = -7j for j = 

1, o-fe = Cfe for fc = 1, . . . ,r; 

(2) Construct bases, V r and W r , that satisfy (4). 

(3) Repeat, while (relative change in {<Ji} > tol) 

(a) A r = W r T AV r and E r = W r T EV r 

(b) Solve the eigenvalue problem A r Xj = AjE r Xj 
and assign Oj < \j for j = 1, . . . , r. 

(c) Update V r so that Range (V r ) = 
span{(<7iE- A) _1 b, ■■• , (oyE - A) _1 b}. 

(4) A r = W r T AV r , E r = W r T EV r , b r = W r T b, and 



these points; instead we use the remaining r degrees of 
freedom to reflect the weight information W(s) and also 
to eliminate terms from the first component of the error 
term. The error expression (9) shows that interpolation 
errors are multiplied by the residues (pi and ipk- Hence, 
we use the remaining r variables to eliminate terms in 
the first and third components of the error expression 
corresponding to the dominant residues (pk and ipk- Note 
that in several cases, such as in the controller reduction 
problem, the state-space dimension of the weight will be 
of the same order as that of G, 0(p) w 0(n). We mea- 
sure dominance in a relative sense; i.e., normalized by 
the largest (in amplitude) <f)k and ipk in every set. More 
details on this selection process can be found in Section 3 
where several examples are used to illustrate these con- 
cepts. Note that one never needs to compute a full eigen- 
value decomposition to obtain the residues of G(s) and 
W(s). Since only a small subset of poles is needed, one 
could use, for example, the dominant pole algorithm pro- 
posed by Rommcs [13] which computes effectively those 
eigenvalues that correspond to the dominant residues 
without requiring a full eigenvalue decomposition. 



Upon convergence of Algorithm 1, cr„- 



-Xj for j 



1, . . . , r; G r interpolates G at these points, and the sec- 
ond sum in (9) is eliminated. W r is unchanged through- 
out, so G r interpolates G at r (aggregated) dominant 
poles of G and W, eliminating v and w terms from the 
first and third sums in (9), respectively. Examples in 
Section 3 illustrate the effectiveness of this approach. 

3 Numerical examples 

We provide two examples related to controller reduction. 
$W and vf/W denote the set of normalized residues of 
G(s) and W(s), respectively. 



3.1 A building model 

The plant, P, is linearized a model for the Los Angeles 
University Hospital, and has order 48; see [4] for details. 
An LQG-based controller, G, of the same order, n = 
48, is designed to dampen oscillations in the impulse 
response. The ten highest normalized residues of G(s) 
and of W(s) are: 

$ (JV) = [ 1.0000 1.0000 0.0286 0.0286 0.0088, 
0.0088 0.0080 0.0080 0.0060 0.0060 ] 

¥ N) = [ 1.0000 1.0000 0.8416 0.8416 0.3935, 
0.3935 0.2646 0.2646 0.0951 0.0951 ] 

There is a significant drop in $( JV ) values after the second 
entry, so we take the first two residues of G as dominant. 
vj/W remains at roughly the same order until the 9 th 
entry. Thus, we choose v = 2; and m = r — v = r — 2 for 
a given reduction order, r. To illustrate the effect of this 
dominant pole selection, we apply W-IRKA, varying v 
from to r. Tables 1 below lists the resulting weighted- 
%2 errors for three cases: r — 12, r — 14, and r = 16. 

r = 12: 



v/vj 


12/0 


10/2 


8/4 


6/6 


4/8 


2/10 


0/12 




1.4021 


1.1433 


0.6548 


0.6863 


0.3576 


0.2181 


0.2853 








r = 


= 14 : 









v /-co 


14/0 


12/2 


10/4 


8/6 


6/8 


4/10 


2/12 


0/14 




1.4734 


1.3436 


0.6477 


0.3019 


0.1538 


0.1425 


0.1351 


0.2224 



16 : 



u/zu 


16/0 


14/2 


12/4 


10/6 


8/8 




1.4206 


1.1934 


0.7258 


0.2898 


0.1917 


u/zu 


6/10 


4/12 


2/14 


0/16 






0.1221 


0.1154 


0.1309 


0.1388 



Table 1 

Weighted-?^ error as v and vj vary 

The weighted- "H 2 error decreases as we take more dom- 
inant poles of W(s) over those of G(s); suggesting the 
importance of the residues of W(s) in the error expres- 
sion (9). Choosing v=2 is the best choice for most cases. 
Tables 1 illustrate that while the weighted error initially 
decreases as v decreases, it starts increasing when v < 2, 
justifying the choice v = 2. For the case of r = 16, sim- 
ilar observations hold Although v = 2 is not the opti- 
mal choice when r — 16, the error for v = 2 is nearly 
smallest, making v = 2 still a very good candidate for 
W-IRKA. These numerical results support the idea of 
choosing v and zu according to the decay of the nor- 
malized residues. Even though this choice seems to yield 
small weighted errors, there may be variations that are 
even better. The residues are multiplied by quantities 
such as W(— Xi)W(Xi), so one might consider incorpo- 
rating these multiplied quantities as well. 

A satisfactory reduced-order controller should not only 
approximate the full-order controller, but also provide 
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the same closed-loop behavior as the original controller. 
Let T and T r denote the full-order and reduced-order 
closed-loop systems, respectively: T corresponds to the 
feedback connection of P with G; and T r to the feed- 
back connection of P with G r . Figure l-(a) depicts the 
amplitude Bode plots of G and G r for r = 14 obtained 
with v = 2. G r is an accurate match to G. Figure l-(b) 
shows that the reduced-closed loop behavior T r almost 
exactly replicates T. 




10 1 10° 10 1 10 2 1 3 1 4 



freq (rad/sec) 



-30 




freq (rad/sec) 



Fig. 1. Bode Plots (a) Full and reduced controller (b) Full 
and reduced closed-loop system 

We now compare W-IRKA with Frequency Weighted 
Balanced Truncation (FWBT) and IRKA of [9] for the 
(unweighted) T~ii problem. Comparison with IRKA is 
included to illustrate the importance of including weight- 
ing in the T^-based model reduction process. We vary 
the reduction order from r — 10 to r — 20 in incre- 
ments of 2, and compute weighted Hoc and %2 errors 
for each case. We use v = 2 for all cases even though 
it might not the best choice for W-IRKA. Results are 
listed in Table 2. Note that for every r value, W-IRKA 
outperforms FWBT with respect to the weighted-?^ 
norm. This might be anticipated since W-IRKA is de- 
signed to reduce the H2 error. But W-IRKA outper- 
forms FWBT with respect to the weighted-'H co norm as 
well in all except the r = 18 case. This is significant since 
balanced truncation approaches generally yield small 
%oo norms. This behavior is similar to the behavior of 
IRKA where one often observes that IRKA consis- 
tently yields satisfactory T-L^ approximants as well [9]. 
Note that for r = 10, the reduced-order controller due 
to FWBT fails to produce a stable closed-loop system. 
Table 2 also illustrates that W-IRKA significantly out- 
performs IRKA in terms of the weighted error norms. 
This is what we have expected since unlike W-IRKA, 
IRKA is tailored towards the unweighted H2 model re- 
duction problem. This becomes clearer after inspecting 
Table 3, which shows that, in terms of the unweighted er- 



ror HG-GVH-H^IRKA outperforms W-IRKA. Thus, 
while G r from IRKA is a better approximation to G in 
an open-loop sense, once the weight is taken into consid- 
eration, W-IRKA does what it is designed for, leading 
to a smaller weighted error. 



r 


10 


12 


14 


16 


18 


20 


FWBT 


1.409 


0.5286 


0.0723 


0.0811 


0.0498 


0.0830 


W-IRKA 


0.9175 


0.1562 


0.0723 


0.0721 


0.0722 


0.0516 


IRKA 


1.4032 


0.4837 


0.1335 


0.0987 


0.1194 


0.1293 



1^ ^ll«oo(W) 



r 


10 


12 


14 


16 


18 


20 


FWBT 


2.1080 


1.1723 


0.1415 


0.1386 


0.1214 


0.1310 


W-IRKA 


0.6677 


0.2180 


0.1351 


0.1309 


0.1028 


0.0956 


IRKA 


1.9540 


0.8620 


0.2102 


0.1286 


0.2066 


0.2317 



\\G L 'r||- H2 ( 1 4/) 

Table 2 

Comparison of W-IRKA, FWBT, and IRKA 



r 


10 


12 


14 


16 


18 


20 


W-IRKA 


0.8152 


0.2750 


0.3679 


0.4078 


0.1274 


0.0518 


IRKA 


0.1062 


0.1168 


0.0478 


0.0513 


0.0123 


0.0082 



IIG-Grllw^/IIGH^ 

Table 3 

Comparison of W-IRKA and IRKA: Unweigthed error 
3.2 International Space Station 12 A Module 

The plant, P, is a model for the International Space 
Station 12A Module with dimension 1412. It is lightly 
damped and its impulse response exhibits long-lasting 
oscillations. A state-feedback, full-order, observer-based 
controller of order n — 1412 is designed to dampen these 
oscillations. The decay rate of the first 50 normalized 
residues $W and are shown in Figure 2- (a). While 
there is almost a two order-of-magnitude drop in $W 
between the third and fourth components, vJfW contin- 
ues to stay significant. Hence, we take v = 3 and re- 
duce order from n = 1412 to r — 60 using W-IRKA. 
For comparison, we also apply FWBT. We denote the 
resulting reduced-order closed-loop systems due to W- 
IRKA and FWBT by T r and T fwb t, respectively. Note 
that Tf w bt was unstable for r — 60. Indeed, r = 88 is the 
smallest order FWBT-derived reduced controller that 
lead to a stable closed-system. All FWBT-derived G r 
are stable; however for r < 88 when G r is connected to 
P, the resulting Tf w bt is unstable. Hence, we compare 
below the r = 60 case for W-IRKA with the r = 88 
case for FWBT. In Figure 2-(b), we plot the absolute 
value of the errors in the impulse responses due to both 
methods. W-IRKA outperforms FWBT even with a 
lower-order controller. We also simulate both T and T r 
for a sinusoidal input of u(t) = cos(2t). Results 2-(c) 
illustrate the superior performance of W-IRKA even 
more clearly. 
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Fig. 2. (a) Decay of the normalized residues (b)-(c) Compar- 
ison of W-IRKA and FWBT using closed-loop responses. 
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5 Conclusions 

We have presented new formulae for the weighted- Tt^ 
inner product and norm that explicitly reveal the contri- 
bution of poles and residues both of the full-order model 
and of the weight. One of the major consequences of 
this new representation are new interpolatory optimal- 
ity conditions for weighted-?^ approximation. Based on 
derived weighted-?^ error expressions, we have intro- 
duced an approach for producing high-fidelity weighted- 
H2 reduced models. The effectiveness of this approach 
has been illustrated with two examples. 
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